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Abstract 



We make an analytic and numerical study of leptogenesis in the framework 
of the (Supersymmetric) Standard Model plus the see-saw mechanism with 
O '• ci U{1) family symmetry and single right-handed neutrino dominance. In 

presenting our analytic and numerical results we make a clear distinction be- 
tween the theoretically clean asymmetry parameter ei and the baryon asym- 
^ ■ metry Yg- In calculating Yb we propose and use a fit to the solutions to 

jL . the Boltzmann equations which gives substantially more reliable results than 

D \ parametrisations previously used in the literature. Our results show that 

"^ ' there is a decoupling between the low energy neutrino observables and the 

leptogenesis predictions, but that nevertheless leptogenesis is capable of re- 
^ ■ solving ambiguities within classes of models which would otherwise lead to 

j^ . similar neutrino observables. For example we show that models where the 

dominant right-handed neutrino is the heaviest are preferred to models where 
it is the lightest and study an explicit example of a unified model of this type. 



1 Introduction 

Leptogenesis is an interesting mechanism which has been proposed to gen- 
erate the observed baryon asymmetry of the Universe (BAU) [|l], |^. The 
mechanism involves the out-of-equihbrium decay of a heavy right handed 
neutrino Nr. The net lepton number L produced in the decay is then re- 
processed into baryon number B by anomalous (B+L) violating sphaleron 
interactions, which otherwise conserve (B-L) 0. 

The advantage of this mechanism is that the same physics that allows the 
right handed neutrinos to decay into light leptons is also responsible for a 
see-saw neutrino mass matrix |Q. This point of view has been strengthened 
by the latest experimental data on the solar neutrino problem by SNO 
and Super-Kamiokande |^ which, when combined, now seems to confirm 
the existence of a solar neutrino mass scale, and suggests active neutrino 
oscillations based on either the LMA or the LOW solution 0. This in turn 
gives impetus to the see-saw mechanism. Combining the see-saw mechanism 
with the experimental data [§, ^ seems to favour scales for right handed 
neutrino masses M^ in the range 10^ - 10^^ GeV. There have been many 
studies of leptogenesis, all based on different models, for example left-right 
symmetry, SO(IO), and so on [^. 

In this paper we study leptogenesis in the framework of the (Supersym- 
metric) Standard Model plus the see-saw mechanism with single right-handed 
neutrino dominance §], [|10[- SRHND is useful for both the LMA and the 



LOW solution |T0[ since it leads to a natural neutrino mass hierarchy in 



the presence of large mixing angles, and gives results which are stable un- 



der radiative corrections |T^. This provides a relatively model independent 
approach which applies to a large class of models with a natural hierarchy 
of neutrino masses [|12|. Indeed in the case of the LOW solution, SRHND 



is almost inevitable in order to maintain the large neutrino mass hierarchy 
present in this case. 

Within the SRHND framework we generalise previously presented an- 
alytic estimates for the mixing angles to the complex domain, and present 
new analytic results for the leptogenesis asymmetry parameter ei and discuss 
the insights which this leads to. We then introduce a f/(l) family symme- 
try |jl3[ and discuss our numerical approach to models of this kind. Our 



analytic results above are supported by the detailed numerical analysis of 
various texture models. Texture models involve unknown coefficients mul- 
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tiplying the expansion parameters, which imphes some level of uncertainty 
in the predictions. In order to quantify this we perform a numerical scan 
over the unknown coefficients, to obtain distributions for predictions of neu- 
trino masses, mixing angles as well as the predictions for ei and the baryon 
asymmetry Yg, for different classes of models. In presenting our analytic 
and numerical results we make a clear distinction between the theoretically 
clean asymmetry parameter ei and the baryon asymmetry Yb- In calculat- 
ing Yb we propose and use a fit to the solutions to the Boltzmann equations 
which gives substantially more reliable results than parametrisations previ- 
ously used in the literature. Using the numerical approach, supported by the 
analytic estimates, we then discuss two important aspects of leptogenesis, 
namely leptogenesis decoupling and leptogenesis discrimination. 

We demonstrate exphcitly that there is a decoupling between leptogene- 
sis and the experimentally measurable neutrino parameters. Although such 
a result may be inferred by comparing the results from different individual 
models which have been proposed in the literature, the present paper repre- 
sents the first attempt to systematically demonstrate this within a framework 
(SRHND) which can be plausibly applied to many different models. To sup- 
port the decoupling claim we present examples of classes of models which 
give the same measurable neutrino parameters but have very different values 
for the CP asymmetry ei. Leptogenesis decoupling implies that there is no 
relation for example between the size of the solar neutrino angle or MNS 
phase and the baryon asymmetry predicted by leptogenesis. 

On the other hand we show that leptogenesis is capable of discriminating 
between different models and thereby resolving ambiguities within classes 
of models giving the same low energy predictions. For example leptogenesis 
may resolve the ambiguity as to whether the dominant right-handed neutrino 
(the one chiefly responsible for the atmospheric neutrino mass in hierarchical 
models) is the heaviest or the lightest of the right-handed neutrinos. We show 
that within a standard hot big bang universe the models where the dominant 
right-handed neutrino is the heaviest are preferred and are more consistent 



with the gravitino constraint on the reheat temperature Tr < 10 GeV W5 . 

In section 2 we introduce our conventions, especially the use of the diag- 
onal charged lepton and right-handed neutrino basis, the see-saw mechanism 
and the MNS matrix in this basis, and the standard model leptogenesis for- 
mulae in this basis. In calculating the baryon asymmetry Yb in section 2.3 
we present and use a new flt formula based on a Boltzmann analysis. In 



section 3 we give our analytic results based on SRHND for the MNS pa- 
rameters and leptogenesis, which give important insights into the numerical 
results which follow. In section 4 we discuss our numerical approach to f/(l) 
family symmetry models. Section 5 is a discussion of the decoupling feature 
of leptogenesis based on the calculation of the asymmetry parameter ei. In 
section 6 we discuss the calculation of Yb for the models where the dominant 
right-handed neutrino is the lightest, and show that such models are not con- 
sistent with a standard hot big bang scenario. In section 7 we then discuss 
models where the dominant right-handed neutrino is the heaviest and show 
that such models can lead to successful leptogenesis. Section 8 concludes the 
paper. 

2 Conventions 

2.1 The Diagonal Charged Lepton and Right-Handed 
Neutrino Basis 

To fix the notation we consider the Yukawa terms with two Higgs doublets 
augmented by 3 right-handed neutrinos, which, ignoring the quarks, are given 

by 



c 



yuk 



+ H.C. 



ij d i j ij u i j ' r\ RR i j 

where e^b = — €&„, ei2 = 1, and the remaining notation is standard except that 
the 3 right-handed neutrinos iV^ have been replaced by their CP conjugates 
iVf and we have introduced a singlet field S whose vacuum expectation value 
(VEV) induces a heavy complex symmetric Majorana matrix Mjm =< S > 
Yrr. When the two Higgs doublets get their VEVS < H^ >= V2, < H} >= 
Vi we find the terms Q 

Cyuk = viY.'jE.E] + v2Y,'<n,n;; + ^mUrN^n;; + h.c. (2) 

Replacing CP conjugate fields we can write in a matrix notation 

Cyuk = ElViY-*Er + NlV2Y'^*Nr + ]^NIM*^^Nr + H.c. (3) 



^In the case of the standard model we must replace one of the two Higgs doublets by 
the charge conjugate of the other, Hd — H^. 



It is convenient to work in the diagonal charged lepton basis 



diag(me,m^,m^) = VeiViY'' Vj^ (4) 



and the diagonal right-handed neutrino basis 



diag(Mi, M2, M3) = V,^M*rrVJ^ (5) 

where VeL-, VeRi Vur are unitary transformations. In this basis the neutrino 
Yukawa couplings are given by 

F- = V,lY''*VJ^ (6) 

and the Lagrangian in this basis is 

'Cyuk = (eL/2LfL)diag(me,m^,m^)(eij/ifir/j)'^ 

+ {l7eLl^^,Ll7rL)Y''v2{NniNR2NR3f 

+ iNR^NR2Nns)dmgiMi, M2, M3)(Ar^iAr^2A^R3)'^ + H.c. (7) 

2.2 The See-Saw Mechanism and the MNS Matrix in 
this Basis 

The light effective left-handed Majorana neutrino mass matrix in the above 
basis is 

ruLL = t^a'^^'^diaglMrS M^\ M^')Y''^ (8) 

Having constructed the complex symmetric light Majorana mass matrix it 
must then be diagonalised by, 

VuLmLLV^L = diag(|mi|, |m2|, jmal) (9) 

where V^l is a unitary transformation and the neutrino mass eigenvalues are 
real and positive. The leptonic analogue of the CKM matrix is the MNS 



matrix defined as |n| 



Umns = KlKl, (10) 

where in the diagonal charged lepton basis VeL will only consist of a diago- 
nal matrix of phases, VeL = Pe corresponding to the charged lepton phase 
freedom, 

/i - Pe /i (11) 



where 

/ e^"^! \ 

Pe = e'*^ (12) 

V e*<^^ j 

These transformations leave the charged lepton masses real and positive, and 
enable three phases to be removed from the unitary matrix V^^, so that Umns 
can be parameterized in terms of three mixing angles Oij and three complex 
phases (5jj, by regarding it as a product of three complex Euler rotations, 



where 



Umns — U23U1SU12 (13) 

/ 1 

U23= C23 8236-^'^^ I (14) 

V -3236'^^' C23 

( Ci3 5136-^^1=' 

t/i3 = 1 I (15) 

V -Si3e*'^i3 Ci3 

/ C12 Si2e-*^i2 
U12 = -sue''^' cu I (16) 

V 1 

where Cij = cos% and Sij = sin 6ij. The resulting MNS matrix is: 

C12C13 Si2Ci3e-*^" sige-**" 

-Sl2C23e'^'^ - Ci2S23Sl3e*^''""*''^ C12C23 - Si2S23Sl3e'^-^''+^''-^''') S23Cl3e-'^' 
5125236*^'^'='+*''^ - Ci2C23Sl3e*'^'" -Ci2S23e'''^=^ - Sl2C23Sl3e'^'''^"*"^ C23C13 

(17) 

The Dirac phase which enters the CP odd part of neutrino oscillation 
probabilities is given by 

5 = 5i3- fe - 5i2. (18) 

2.3 Leptogenesis in this Basis 

CP violation in the decay of the lightest right-handed neutrino Njn comes 
from the interference between the tree-level and one-loop amplitudes B & 



p!6| , |T^. The CP asymmetries given by the interference with the one-loop 
vertex amphtude are in the SM [0, ||: 



ei 



T{N, 



Rl 



L, + H2) 



r(ivii 



Lj^ + Hi 



T{Nbx 
1 



.Ml 



87r(rjF, 



Y.Im( iXX 



1^)11 i^i 



u)li 



.Ml 



where 



/(^) 



'1 + x) In 



1 +x 



X 



^<M?' +»(«?' 



g[.x) 



X 



X 



(19) 



(20) 



where f{x) arises from the interference between the tree level decay and the 
vertex correction, while g{x) is due to the interference with the absorptive 
part of the one-loop self-energy, which can in principle be much larger if 
the right-handed neutrinos are almost degenerate |I6| , |T^. Assuming that 
Ml <^ M2 <^ M3, we have approximately ||T8[| , 



ei 



167r(rji;)n t^^ 



Y.Imi\{YX)i 



Ml 



(21) 



In the Supersymmetric SM the result for ei is twice as large as in Eq.^ due 
to the extra SUSY degrees of freedom in the diagrams. 

The lepton asymmetry Yl of the universe created by this mechanism can 
be written as 

Yl = d^ (22) 
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where ei has been defined above, g* counts the effective number of degrees 
of freedom, for the SM g* = 106.75 while for the Supersymmetric SM g* = 
228.75 [^ and d is the dilution factor which takes into account the washout 
effect produced by inverse decay and lepton number violating scattering. To 
calculate d one has to solve, in principle, the full Boltzman equations, which 
can be done numerically 0, |18| . 



However, many authors, for examples see ^^, have used simple approxi- 
mated solutions to the Boltzman equations expressed as 



d 



0.24 
k{lnkf/^> 



10 < A; < 10^ 



(23) 



d = — 
2k 

d=l 



l<k<10 
0< k<l 



(24) 
(25) 



Recently, Nielson and Takanishi ||2T| suggested a slight modification of 
eqs (^3])-(^), namely, 

0.3 



d 



d 



k{lnkfl^ 
1 



2VF+9 
Here the parameter k is given by. 



10 < fc < 10^ 
A; < 10 



k 



Mr 



{Yiyu] 
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1.7 X Sttv^ Ml 



(26) 
(27) 

(28) 



where Mp is the Planck mass. Physically A; ~ 1 represents the desirable re- 
gion where the couplings of the right-handed neutrinos are sufficiently strong 
for them to be copiously produced from particles in the thermal bath, but 
sufficiently weak for them to decay satisfying the out-of-equilibrium condi- 
tion. 

We have compared eqs (p6|)-(p7D to the full numerical solution to the 
Boltzman equations as plotted in fig. 6 of |19|. Buchmiiller and Pliimacher 
[ p!9|] use as a variable 



mi 



: {YX 



u II 



Ml 



1.1 X 10"'^ A; eV 



(29) 



The result of this comparison is shown in Fig. |1]. As one can see, the 
approximate formulas, eqs (|26|)-(p7D, are a valid approximation for values of 
fc ~ 1 and for small values of Mi < 10^° GeV. For smaller values of fc <C 1 
or fhi <^ 10"'^ eV the approximation formulas are clearly not valid since 
they do not take into account the production suppression apparent in the 
full treatment using the Boltzmann equations. For larger values of fc 2> 1 or 
fhi 3> 10~^ eV (and larger values of Mi) the approximation formulas are also 
clearly not valid since they do not take into account the steep suppression due 
to the out-of-equilibrium condition being violated which is again apparent in 
the full treatment using the Boltzmann equations. From Fig. ^ it is obvious 



that in this case the analytic approximation seriously underestimates the 
suppression of d by orders of magnitude. 

For this reason we have devised a purely empirical fit formula for the 
exact solution to the Boltzman equations, which can be written as, 



(a) Logi^ids^L) = 0.8 * Logw{m^) + 1.7 + 0.05 * Logi^{Ml^) (30) 
(6) Logi^{dB-L) = -1.2- 0.05 *Lo(7io(Mi^°) (31) 

{c) Logi^{dB-L) = -(3.8 + Lo(?io(Mii°))*(Lo(7io(mi)+2) (32) 

- (5.4 - ^ * Lo(7io(Mi))=^ - ^ 

where MY' = Mi/lO^'^GeV, fhi is in units of [eV]. In implementing this fit it 
is always the smallest of (a)-(c) which is taken. The results of this fit are 
superimposed onto the exact curves taken from |]19| in Fig. |l]. 

Obviously eqs. (^)- ( ^2|) reproduce the exact results considerably better 
than eqs (|26| ) - (|27|) . However, in |T^ the authors assumed that right-handed 



neutrinos are hierarchical. Also we have to assume that for small values of 
ei the dilution function does not depend on ei. Thus, eqs (pOl)- (|3^) are still 
approximate, and may not be valid if these conditions are violated. Moreover, 
for large values of Mi, say Mi ~ 10^^, we have to rely on an extrapolation 
beyond values of Mi used in the fit. 

Note also that since in the SUSY SM both ei and g* are twice as large 
as in the SM, the two effects tend to cancel in the estimate of Y^. Also the 
approximations for d over the above range of rhi are valid for either the SM or 
the SUSY SM 12^]. Therefore the results we present are approximately valid 



for either the SM or the SUSY SM, although for definiteness we consider the 
SM from now on. 



Due to sphaleron effects Yl finally is related to Yb approximately via p2 

Yb = ^—Yl (33) 

a — i 

where 



•^We fit ds-L since authors of [Q plot Yb-l, where ds-L is related to d via d 
(1 — a)dB^L, where a is defined in Eq.|34 
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Figure 1: The logarithm of the dilution function (Ib-l versus the logarithm of 
rhi for different values of the lightest right-handed neutrino mass Mi = 10^, 
10^°, 10^2 Qgv (from right to left). The thick line is the solution to Eqs.(^)- 
(pTl). The dashed, medium full, dotted curves represent cIb-l for Mi = 10^, 
10^°, 10^^ GeV extracted from Fig. (6) of |TP|, based on an exact numerical 
solution of the Boltzman equation. The thin solid curves with plateau regions 
are the approximately fitted values of ds-L following Eqs. 
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Figure 2: The approximate solutions of the fit function for (Ib-l from 
Eqs.(§D])-(§g)) for (from right to left) Mi = 10^ 10i\ 10^^ 10^^ GeV. For 
further discussion see caption in Fig.|l[ 



Here Np is the number of families and Nh the number of Higgs doublets. 
In the SM a ~ 1/3. Experimentally Yb is expected to be in the range 
Yb = (ub - nB)/s ~ (0.5 - 1) X lO-^" @, Q. 



3 Analytic Estimates 

3.1 MNS Parameters from SRHND 

In the basis used in this paper where the charged leptons are diagonal, and the 
right-handed neutrinos are diagonal, we write the neutrino Yukawa matrix 

as 

(a' a d \ 

h' h e \ (35) 

d c f I 

where the LR notation means that the second and third columns of Y^ cor- 
respond to the second and third right-handed neutrinos. We use the phase 
freedom of the charged lepton masses in Eq.[Tl| to make the couplings to the 
third right-handed neutrino d, e, / real and positive, leaving a, 6, c, a', 6', c' 
complex. 

We write the diagonal (real, positive) Majorana masses in this basis as 



M, 



RR 



( X' 











X 





I 





Y 



(36) 



Then using the see-saw formula for the light effective Majorana mass matrix 
rriiL = fly^M^^y "^ (valid for complex couplings) we find the symmetric 



matrix, 



rriLL 



I oP_ _|_ of. _|_ c^ g'b' _|_ o^ _|_ rfe o-' (^ i gc i d/ \ 

' X' ^ X ^ Y X' ^ X^ Y X' ^ X ^ Y > 

fe^ifefief, h'c' . be . ef 

X' ~^ X ~^ Y X' ~^ X ~^ Y 

V ■ ■ x^ + X + T^ / 



(37) 



So far the discussion is completely general. In order to account for the 
atmospheric and solar neutrino data many models have been proposed |]12 



based on the see-saw mechanism [Q. One question which is common to all 
these models is how to arrange for a large mixing angle involving the second 
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and third generation of neutrinos, without destroying the hierarchy of mass 
sphttings necessary to account for the solar and atmospheric data. Assuming 
6*23 ~ 7r/4 one might expect two similar eigenvalues m2 ~ m^, and then a 
hierarchy of neutrino masses seems rather unnatural. 

For our analytic estimates, we assume for simplicity that the first right- 
handed neutrino X' contributions are insignificant compared to the second 
right-handed neutrino X contributions, 

« T^ (38) 



X' X 

Then one way to achieve a natural hierarchy is to suppose that the third 
right-handed neutrino contributions are much greater than the second right- 
handed neutrino contributions in the 23 block of uill 10 , 



Y X 

This implies an approximately vanishing 23 sub determinant, 

c^etK,]23=(- + ^)(^ + -)-(^ + ^) ^0 (40) 



and hence 



1 1 


(ef be 


Y x)[y x) 


[y ^ X 


1^12/1713, < 1 





(41) 



Thus the assumption in Eq.^ that the right-handed neutrino Y gives the 
dominant contribution to the 23 block of rriLL naturally leads to a neutrino 
mass hierarchy. This mechanism is called single right-handed neutrino dom- 
inance (SRHND)[0. In the limit that only a single right handed neutrino 
contributes the determinant clearly exactly vanishes and we have 1712 = ex- 
actly. However the sub-dominant contributions from the right-handed neu- 
trino X will give a small finite mass 1712 7^ as required by the MSW solution 
to the solar neutrino problem. 

Assuming SRHND as discussed above, we may obtain a simple estimate 
for the third neutrino mass: 



2id' + e' + f 



^3 - v^^ '- (42) 
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Note that mi^2 are determined by parameters associated with the subdomi- 
nant right-handed neutrinos and so are naturally smaller. Given the SRHND 
assumption in Eq.|3^ we see that we have generated a hierarchical spectrum 

In order to obtain the MNS parameters we must diagonalise rriLL as in 
Eq-i 

VyLmLiV^L = diag(|mi|, jma], Img]) (43) 

where we write V^l as a product of complex Euler rotations of the form of 
Eqs.|13|,|l^,|T5|,|16|, together with diagonal phase matrix 



/ e^^i/2 \ 

e^'^2/2 
y e^^3/2 



(44) 



which is required to remove the phases in rrii = \mi\e^^% 

y.L = Plululul (45) 

Thus VuL contains the 3 angles and 6 phases of a general unitary matrix. 
However in the basis where we have chosen the couplings c?, e, / to be real, 
ms is given in Eq.^ and 03 is zero to leading order. 

In order to bring the MNS matrix into the form in Eq.^, additional 
charged lepton phase rotations are required as in Eqs.|TD|, so that we have 
finally 

UmNS = PeU23Ul3Ul2Pu (46) 

where Pe is a diagonal matrix of phases as in Eq.|T2|. Note that the angles 
involved in the Uij are the same as those in the Uij in Eq.|13], 6ij = 6ij, but 
the phases will be different, 6ij ^ 6ij, due to the non-zero phases in Pe,Pu- 

Since the couplings d, e, / are real, we find that the previous estimates 
based on SRHND are still vahd [0 

tan^23~4 (47) 
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where the associated phases are approximately zero 

.^23 ^ ^13 ^ (49) 

By a suitable choice of parameters e = / ^ rf it is possible to have maxi- 
mal ^23 suitable for atmospheric oscillations, while maintaining a small ^13 
consistent with the CHOOZ constraint ||23|| . 

To determine U12 is quite complicated in general, but in the physically 
interesting cases where 612 is near maximal 612 ~ 7r/4 we find the simple 
analytical results 

tan0i2^y2^:^ (50) 

\b — c\ 

^12 ~ (Pb-c - (pa (51) 

where 06_c = arg{b — c) and (pa = arg{a). In the simple example that the 
phases in P^, Pi, are zero, the observable Dirac phase in Eq.[l^ is given in 
Eq.^. In general the Dirac phase will involve a more complicated combina- 
tion of phases. 

3.2 Leptogenesis in SRHND 

In leptogenesis it is generally the lightest right-handed neutrino which decays 
to produce lepton number, where we use the notation that Mi is the lightest 
right-handed neutrino, M3 is the heaviest right-handed neutrino and we as- 
sume Ml -C M2 <C M3. In the notation of the previous subsection where Y 
is the dominant right-handed neutrino there are two physically distinct cases 
to consider: 

(a) F < X < X' (i.e. Y = Mi, X = M2, X' = M3) 

(b) X' < X < F (i.e. X' = Mi,X = M2, Y = M3) 

In other words the dominant right-handed neutrino may either be (a) the 
lightest, or (b) the heaviest right-handed neutrino, and both cases must be 
considered. 

It is also worth emphasising that there is no generation ordering implied 
by the results in the previous subsection (or those in @,ll3). In other words 



the dominant right-handed neutrino Y may be associated with the third, 
second or first generation, by a simple reordering of the columns of Yi,. Due 
to the hierarchy of charged lepton masses, it is meaningful to associate the 
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first row of Yi, with tlie first generation, the second row of Y^, with the sec- 
ond generation, and the third row of Y^, with the third generation. However 
the physical neutrino mass matrix rriiL is invariant under the operation of 
exchanging the columns of Y^, along with the ordering of the right-handed 
neutrinos in Mjir, so the SRHND results apply quite generally to all gen- 
eration orderings of the right-handed neutrinos p, [0. Physically if the 
Yukawa couplings e, / are of order unity, then it may be natural to associate 
Y with the third generation. However if the couplings e, / ^ 1 then it may 
be more natural to associate Y with the second generation, and re-order the 
matrices by interchanging of the second and third right-handed neutrinos in 
Y^ and Mrr. 

Returning to the leptogenesis asymmetry parameter in Eq.|2T|, for case 
(a), where the dominant right-handed neutrino mass Y is the lightest, using 
the SRHND results of the previous subsection, we find 

'^''^~ ^ ^^\in(205+c)|& + c|2 (52) 



327r VX 

while for case (b), where the dominant right-handed neutrino mass Y is the 
heaviest, we find 



^) 




\b' + c' ^^ 

^'|2 + |6'|2+ |c'l2 



sm{2<Pb'+c')e , ,|n , |^,|n , , ,|n (53) 



where (ph+c = arg{b + c) and (pb'+c' = arg{b' + c'), and we have used the fact 
that 7712 "^ '^s in obtaining Eq.|53|. Q 

Are these values of ei of the correct order of magnitude? We may use 
m^ ~ 5 X 10~^eV and m^ ^ ^^2 y^ i^ Eq.^ and the crude order of magnitude 
approximation for 1712 ~ |6 — cpf|/M2, to obtain 

.« ~ si„(2,.«)10- (;^) (^^) (54) 



ef> ~ siu(2^.,,)10-= (^ji^) (55) 



■'it is also apparent that the phases which are relevant for leptogenesis in both cases 
are not identical to the Dirac phase which even in the simple example that the phases in 
Pe,Pu are zero, is given as in Eq.Qq as 5 ~ (pa — 4'b-c- In general the Dirac phase will 
involve a more complicated combination of phases still. 
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The results in Eqs.|5^, p5| express ei in terms of the hghtest right-handed 
neutrino mass in each case. Since e^ is suppressed relative to e^ by a factor 
of 1112/111^ (which should be 7712/ m^ < 0.1), this implies that the lightest 
right-handed neutrino mass must be at least an order of magnitude larger in 
case (a) than in case (b). 

To understand which of the two cases (a) or (b) is more promising from 
the point of view of leptogenesis it is important to estimate the parameter 
rhi in Eq.^ which controls the dilution factor as shown in Figs. ^^. From 
Eqs.||,|||, 

m'f^ ^ t^ r'^' ' ^\l ' '" ^ (56) 



2i\d' 


+ 


|e2 


+ 


\m 


2 

1 /|2 


+ 


Y 

6f 


+ 


cf) 



mf ^ v t^-^ ' '^'^ ' '-'^ (57) 

In case (a), where the dominant right-handed neutrino is the lightest one, 
the parameter fhi in Eq.^ is approximately equal to the physical mass of 
the heaviest neutrino in Eq.^ which is measured by Super-Kamiokande. 
Thus for these models m^ ~ m^ ~ 5 x 10"^ eV which is generally beyond 
the plateau regions in Figs.|l|, |^, and this leads to the requirement that Y ~ 
Ml < 10^ GeV and a strong dilution suppression d <C 1. However, according 
to Eq.|4|, y ~ Ml < 10^ GeV leads to values of ej"^ < 10"^ which, when 
combined with the dilution suppression (i -C 1, implies from Eq.^ Yb ^ 
10"^'' well below the observed value. 

In case (b) , on the other hand, where the dominant right-handed neutrino 
is the heaviest one, there is no association of fhi in Eq.^ with a physical 
neutrino mass and so this parameter may in principle take smaller values 
closer to the plateau regions, leading to only a mild dilution suppression 
d < 0.1 for a range of lightest right-handed neutrino mass X'. Furthermore, 
as we already remarked, by comparing Eqs.p^, E^ we see that value of ei 



is larger by an order of magnitude than ei . For example if we choose 
X' ~ Ml < 10^ GeV, consistent with the gravitino constraint on the reheating 



temperature Tr < 10 GeV ||I5[, we find ei < 10 which, assuming a mild 



(b) 



^Note that since the dominant right-handed neutrino mass is given by Y ^ e^5.10^'' 
GeV, case (a) requires e ^ 1, whereas for case (b) it is consistent with e ~ 1 providing 
there is a sufficiently large hierarchy in the right-handed neutrino sector. This means that 
in case (a) the dominant right-handed neutrino cannot be associated with the third family, 
whereas in case (b) it may be. 
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dilution suppression d < 0.1, implies from Eq.^ Yb ~ 10^^'^ which is just 
about acceptable. 

We shall later present specific examples with detailed numerical results 
which support the conclusion that leptogenesis prefers case (b) where the 
dominant right-handed neutrino is the heaviest one, at least according to the 
standard hot big bang picture, ignoring effects of infiation. 

4 Numerical Approach to U{1) Family Sym- 
metry Models 



Our numerical results are based on the SRHND models P],[T^, with a U{1) 



family symmetry. The idea of such a symmetry is that the three families 
of leptons are assigned different f/(l) charges, and these different charges 
then control the degree of suppression of the operators responsible for the 
Yukawa couplings, leading to Yukawa matrices with a hiearchy of entries. 



and approximate "texture" zeroes |]T3[. As usual it is assumed that the U{1) 



is slightly broken by the VEVs of some fields 6, 6 which are singlets under the 
standard model gauge group, but which have vector-like charges ±1 under 
the U{1) fiavour symmetry. The U{1) breaking scale is set hy < 9 >=< 9 >. 
Additional exotic vector matter with mass My allows an expansion parameter 
A to be generated by a Froggatt-Nielsen mechanism [^, 

^ = ^ = ^-0-22 (58) 

where the numerical value of A is motivated by the size of the Cabibbo angle. 
Small Yukawa couplings are generated effectively from higher dimension non- 
renormalisable operators corresponding to insertions of 9 and 9 fields and 
hence to powers of the expansion parameter in Eq.^ The number of powers 
of the expansion parameter is controlled by the U{1) charge of the particular 
operator. The lepton doublets, neutrino singlets, Higgs doublet and Higgs 
singlet relevant to the construction of neutrino mass matrices are assigned 
t/(l) charges k, Up, hu = and a. From this starting point one may then 
generate the neutrino Yukawa matrices as in 0. The neutrino Dirac Yukawa 
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matrix is 

/ anAl'i+"il ai2Al'i+"2l ai3Al'i+"='l 

f = a2iAl'2+"il a22Al'2+"2l a23Al'2+"3l | (59) 

V a3iAl'3+"il a32Al'3+«2l a33Al'3+"3l 

where Eq.|5^ may be identified with Eq.^. The heavy Majorana matrix is 

AiiAl2"l+'^l Ai2Al"i+"2+a| ^^^^|„i+„3+<x| 
Yrr = I Ai2Al"2+"i+'^l A22Al2"2+'^l A23Al"2+"3+<x| I (^QO) 

^^3AI"3+ni+a| ^^_^^|„3+„2+a| ^33Al2"3+a| 

where Ajj and Ojj are undetermined coefficients, and A is an independent 
expansion parameter relevant for the right-handed neutrino sector. Q 

The neutrino Yukawa matrices are generated in a particular basis defined 
by the f/(l) family symmetry. This corresponds to the starting basis defined 
by tildes in section 2, and numerically we follow the procedure to go to the 
diagonal right-handed neutrino basis, as outlined there. Note that we assume 
as an approximation that the charged lepton matrix is diagonal with positive 
eigenvalues in the starting basis. In practice this may be approximately 
achieved by a suitable choice of right-handed lepton f/(l) family charges, as 
discussed elsewhere [^,[llO|. 

In our numerical analysis we take account of the fact that the theory does 
not determine the complex coefficients Aij and aij which one has to choose in 
some range. This is not a special feature of the SRHND models, which we are 
focussing on in this paper, but a limitation of texture models based on a U{1) 
family symmetry. Usually one simply assumes that the unknown coefficients 
are of order 0{1) and, therefore, the structure in the Yukawa matrices is given 
by the expansion parameter rather than the coefficients. Our approach to this 
problem is to scan over the unknown coefficients randomly and to construct 
distributions for the various observables of interest. This way we are able to 
determine distributions for masses and mixings of a given model. Given the 
statistical nature of our approach, one question comes immediately to mind: 
What is the correct range of values one should choose for the coefficients? 
Lacking any theoretical background we have chosen for the coefficients the 
interval 



^We are grateful to G.Ross for emphasising that the right-handed neutrino sector is 
controUed by an independent expansion parameter. 
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(^ij J 



Aij =^ [^/X, 1/V\] X 6*"^'^ (pij => [0, 27r] (61) 



It should be noted, that this choice is the minimum requirement for texture 
models to be sensible, simply because any larger variation in the coefficients 
would destroy the texture one originally assumed to be the dominant feature 
of the mass matrices of interest. 

A word of caution might be in order. Obviously the distributions which 
we calculate depend on our choice for aij,Aij. Lacking further theoretical 
support for our choice, we can not evaluate the success of a given model in 
terms of confidence intervals. Instead our method is more minimalistic. We 
will consider a model to be a "good" model, if the main body of the distri- 
bution in a given observable coincides with or is close to the experimentally 
preferred value. Clearly, a model which fails even our simplistic test will fail 
even more badly under a more sophisticated numerical analysis. We would 
like to stress, however, that although the width of the peaks and the de- 
tailed shape of the distributions change under a change of the range of the 
coefficients, the position of the peaks remains nearly invariant. 

In order to be able to compute the expectations for the leptogenesis "ob- 
servable" ei in the different models, our current computation goes beyond the 
one we discussed in a previous paper |]2^ in allowing the coefficients aij, Aij 
to be complex. Since we do not have a theory of phases, we decided to choose 
the (pij in the full interval [0, 2n]. In other words, since we do not know about 
any mechanism suppressing phases effectively in the Yukawa couplings, we 
simply expect that all phases should be large. 

So, our numerical procedure may be summarised as follows. First select 
a particular flavour model defined by a choice of f/(l) charges. Second select 
randomly a set of complex coefficients aij and Aij. Third diagonalise the 
right-handed neutrino mass matrix to yield positive eigenvalues, and express 
the Dirac Yukawa matrix in this basis, as discussed in section 2.1. Fourth 
calculate the see-saw matrix miL and hence the physical neutrino masses and 
the MNS angles and three phases as discussed in section 2.2. Fifth calculate 
the leptogenesis parameters ei and Yb as discussed in section 2.3. Then the 
whole procedure is repeated for a different set of randomly chosen complex 
coefficients aij and Aij, and the results are binned to build up distributions 
of the observable quantities. In the figures we show in the following we use 
random sets of 10*^ matrices for each of the distributions shown. Finally a 



different model is selected corresponding to a different set of U{1) charges and 
the whole procedure is repeated for the new model. We disregard the effect 
of renormalisation group radiative corrections in going from high energy to 
low energy, which has been demonstrated to be of the order of a few per cent 



for SRHND models 11 



5 Leptogenesis Decoupling 

In this section we will discuss leptogenesis decoupling, namely, the fact that 
the leptogenesis observable ei can take any value independent of the low 
energy observables, i.e. masses and mixings. Unfortunately this means that 
measurements of the solar angle or the MNS phase for example does not 
tell us anything about leptogenesis. On the other hand the results in this 
section also demonstrate another aspect of leptogenesis, namely that it can 
be used to resolve the ambiguity between different models which all lead to 
very similar predictions for low energy neutrino observables. In this way 
leptogenesis provides information about the high energy theory which would 
be impossible to determine by the measurement of low energy observables 
alone. 

We will defer the discussion of Yg until the next section and concentrate 
here only on the calculation of ei, since the conversion of the CP asymme- 
try parameter to Yb depends highly on the assumed thermal history of the 
universe, whereas the calculation of ei is theoretically clean. 

In Table |I| we give four examples of models based on different choices 
of flavour charges. For simplicity, we start by assuming that the expansion 
parameter in the right-handed neutrino sector is equal to the Wolfenstein 
parameter A = A, as was assumed in |T0[. Model FCl was discussed analyt- 



ically in |jTO|, where it is seen that it yields a heavy Majorana matrix with 
an off-diagonal structure in the U{1) charge basis. It satisfies the SRHND 
conditions, and has a r^ b,c and so leads to the LMA solution. FC2 also 
has an off-diagonal heavy Majorana matrix, but has a <^ b,c and so leads 
to the SMA solution. |^ FC3 is also taken from [|1^], and is an example 
of a model with an approximately diagonal heavy Majorana matrix in the 
U{1) charge basis. Using the analytic results in |T0[ we find the approxi- 

^The latest data from the SNO collaboration ||] rather strongly disfavours the SMA 
solution 0. 
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Table 1: Flavour charges (FC) for four models, as discussed in the text, and 
approximate expectations for ^23; ^13, ^12 and R for the four different models. 



mate expectations for the experimentally accesible quantities (6*23, ^13, 612 
and R = |Am2i|/|A??i32|) where Am?- = mf — m'j as given in Table |I[ Thus 
FCl is suitable for the LMA solution, FC2 for the SMA solution, FC3 for 
the LMA but with a larger CHOOZ angle than FCl, and FC4 is a model 
without SRHND which is consequently expected to give a larger value of R 
than models FC1-FC3 which all have SRHND. f\ 

Figure 1 shows the distributions for the solar {sq = 4 sin ^^2(1 ~ sin ^^2))) 



atmospheric {sAtm = 4sin6'|3(l — sin6'23)) and CHOOZ {sc = 4sin6'?,(l 



'13 >> 



sin ^^3)) angles as well as for R = \A'm'2i\/\Aml2\ for the four models given 
in table 1. As discussed above, the detailed shape of the distributions is 
different to the one we calculated previously [0] using real coefficients. The 



positions of the peaks of the various distributions, however, did not change 
allowing for complex phases. 

Figure ^ shows the distributions in the leptogenesis observable e for the 
models FC1-FC4. From the figures one might be tempted to think, that 
different low energy observables lead to different values of ei and so might 
be distinguished. This is not true, and we now show that any of the models 
can be modified to give any desired value of ei, while keeping the low energy 
observables approximately unchanged. 

Let us consider as an example the model FC3 discussed above, which 



^As a side remark we mention that for neutrinoless double beta decay, in flavour models 
which make use of the seesaw mechanism, one never expects that the effective Majorana 
neutrino mass {{rriiy) — (toll)ii) measured in double beta decay is exactly zero. However, 
these models produce a hierarchical spectrum of left-handed neutrinos and thus one expects 
{rriiy) to be small. In the models we have studied in this paper, one typically gets (m^) ~ 
10^'^ eV, albeit depending on the model and with a rather larger uncertainty. 
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Figure 3: Theoretical distributions for the predictions of neutrino mass and 
mixing parameters for four selected see-saw models: FCl (full), FC2 (dot- 
dashes), FC3 (thick dots), FC4 (dashes). Matrix coefficients are randomly 
chosen in the interval [\/2A, l/v2A] x e*''^, with =^ [0, 27r]. The vertical axis 
in each panel (deliberately not labelled) represents the logarithmically binned 
distributions with correct relative normalisation for each model, with heights 
plotted on a linear scale in arbitrary units. 
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Figure 4: Plots of e for the four different models of table 1. Plot style for 
different models follows fig. ^ 



21 



predicts the LMA solution and a relatively large CHOOZ angle. Model FC3 
gives (neglecting the coefficients and assuming A = A) the following Dirac 
and Majorana mass matrices: 

(yl/2 ^ y^/2 \ 
A3/2 A Ai/2 (62) 

A3/2 A Ai/2 j 

( 1 Ai/2 A \ 
YK^' = ■ a A3/2 (63) 

which after see-sawing give the following leading order structure for mLL. 

(64) 




Note, that FC3 has a right-handed neutrino mass matrix which is diagonal 
to leading order and it is the lightest (third) right-handed neutrino which 
gives the dominant contribution to niLL. The estimate for the asymmetry 
parameter is given in Eq.|5^, where it is clear that e ~ A^^^. In order to change 
ei we must reduce e. This may be achieved by adjusting the U charges in 
such a way that the Dirac neutrino matrix just gets multiplied by an overall 



scaling factor compared to eq. ^ while the heavy Majorana Yukawa matrix 
remains unchanged. The rescaling of the Dirac Yukawa matrix implies that 
the coupling e is made smaller, and hence the scale of right-handed neutrino 
masses must be reduced in order to maintain the same value of m^. This 
will lead to a different value of ei without changing the other low energy 
observables at all. 

This qualitative conclusion is supported by our numerical results. In table 
y we give sets of charges for variants of the model FC3 of table |l|, which lead 
to a simple rescaling of Y^,, 

YFC3{a,b,c,d,e) ^ Q(l,2,3,4,5)\yFC3 /gg-j 

and hence the scale of right-handed neutrino masses as shown in fig |. 

All of these models were constructed to preserve the low-energy phe- 
nomenology, and in fact we have checked that they lead to identical predic- 
tions for SAtm, Sq, sc and R as FC3. Figure |^ shows the resulting values 
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Table 2: "Variants" of the flavour model FC3 of table |^. All of these models 
give exactly the same distributions for the low energy neutrino observables. 
They differ, however, in their predicted values for the leptogenesis obervable 



e. 



of Ml and e. Note that the lightest right-handed neutrino mass for FC3a, 
FC3b and FC3c in Fig.^ is above the reheat temperature allowed by the 
gravitino constraint [0. This figure explicitly demonstrates that it is possi- 
ble to completely decouple the predictions for leptogenesis from low energy 
observables. 

How well do the analytic estimates for e discussed previously agree with 
the numerical results? In terms of our small expansion parameter A ~ 0.22 
and inserting the flavour charges for the models FC3 (FC3a, FC3b, FC3c, 



FC3d and FC3e) into eq. ^ one finds 



M 



327r 



A3 (A^A^A^Al^Al3) 



(66) 



numerically 3 x 10"^ (2 x 10-^ 7 x 10"^ 4 x 10-^ 2 x lO^^ and 8 x lO"") 
which coincides approximately with the peaks of the distributions in e shown 
in fig. 1^. Recall that model FC3 predicts a right-handed neutrino mass 
matrix with the dominant neutrino being the lightest one (case a, discussed 
in section 3.2). 

Model FC3 produces predictions for low-energy neutrino phenomenology 
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Figure 5: Plots of the mass of the hghtest right-handed neutrino Mi (left) 
and ei (right) for the five different models of table 0. From right to left: 
FC3a-FC3e. 

consistent with the large angle MSW solution of the solar neutrino problem. 
It is interesting to ask whether this solution is the only one for which one 
can decouple ei from the low-energy observables. 

In order to investigate this problem we have constructed variants of FC2, 
predicting a small angle MSW solution to the solar neutrino problem. The 
corresponding charges are given in Table. |^. All models in this table produce 
exactly the same distributions in R and SAtm as model FC2, but lead to 
different predictions for s©, sc and ei as is demonstrated in Fig. ^ Note 
that we have multiplied the distributions for FC2a and FC2b by a factor 
of 1.1, since otherwise the curves would completely overlap in some of the 
variables. 

As can be seen from Fig. ^ models FC2 and FC2a give the same predic- 
tions for Sq and sc, but differ in their predictions for ei. FC2b and FC2c, on 
the other hand, give expectations for Sq and sc which are smaller than the 
one for FC2 by about 1.5 orders of magnitude. Nevertheless, FC2b yields 
values of e which are very similiar to those of FC2. Also FC2c and FC2a 
have very similar expectations for leptogenesis while differing in Sq and sc- 

It is obviously easy to find models differing in their predictions for lepto- 
genesis and at the same time being consistent with SMA MSW. Moreover, 
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Table 3: "Variants" of the flavour model FC2 of table |1|. FC2 and FC2a give 
the same distributions for the low energy observables, with an expectation 
for the solar and CHOOZ angle of order A^, FC2b and FC2c, on the other 
hand, lead to an expectation for solar and CHOOZ angle of order A^. See 

fig. I 

neither the size of the solar nor the size of the CHOOZ angles tell us anything 
about whether leptogenesis is possible or not. 

Finally we have investigated the question whether a special value of R 
determines ei. All the models discussed so far prefer values oi R> 10~^. The 
following assignment of charges defines a model (FC5), which prefers larger 
hierarchies, see fig. ^ 



(/i, k, h, ni, 722, ng, a) = (3, -3, -3, 0, -1/2, 1, 1), 



(67) 



while still keeping the atmospheric and solar angles large (and sc <^ 1)- 
FC5 therefore is consistent with the LOW solution of the solar neutrino 
problem. Nevertheless, as fig. demonstrates FC5 leads to a very similar 
expectation for e as the model FC3b discussed above, which prefers R in the 
range R ~ lO'^ - IQ-^. 

Obviously, any value of R can produce approximately the same order-of- 
magnitude values of ei. 

As a summary it can be stated that there is a decoupling between low 
energy neutrino observables and leptogenesis. We have demonstrated this 
point by constructing a number of different flavour models, which give the 
same predictions for neutrino masses and mixings while differing by huge 
factors in their expectations for leptogenesis. 

On the other hand we have seen that leptogenesis is in principle able to 
resolve the ambiguity between different models which would lead to the same 



25 





10"^ lO-"* 10-3 10-2 10-^ 10° 

sc 



10-' 10-" 10-=^ 10-^ 10- 



Figure 6: Plots of Sq (top left), sc (top right), e (bottom left) and Yb (bottom 
right) for the four different models of table ^. The full line is FC2, the dashed 
line FC2a, the dotted line FC2b and the dash-dotted line FC2c. Note, that 
the distributions for FC2a and FC2b have been multiplied by a factor of 1.1, 
see text. 
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Figure 7: Plots of R (to the left) and e (to the right) for the 2 different 
models FC3b (full line) and FC5 (dashed line). 
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low energy neutrino observables, and which otherwise would be indistinguish- 
able. Therefore leptogenesis is able to provide information about the high 
energy theory which could not be obtained by low energy measurements. 

6 From ei to Yg when the dominant right- 
handed neutrino is the lightest 

While the calculation of ei is straightforward, once a particular model has 
been specified, the calculation of Yb depends crucially on a number of as- 
sumptions made about the early universe. In the following calculation we will 
assume a standard hot big bang scenario in which the maximum temperature 
is higher than the largest right-handed neutrino mass we consider, such that 
the right-handed neutrinos can be thermally produced. This assumption is 
necessary if one wants to employ one of the parameterizations to the full 
solution of the Boltzman equations, see eqs 2^27 and 30-3^, discussed in 
section 2.3. 

Obviously the following discussion will not be valid, if the universe under- 
went a period of inflation with a rather low reheat temperature, as required 
by the gravitino problem. 

Let us first discuss the different values for Yb one obtains using either eqs 
2B| - P7| or our parameterization eqs P(]| - P^ As an example we will concentrate 



on the variants of the model FC3, discussed in the last section. 

Fig. H shows calculated values of Yb using the two different approxima- 
tions. Obviously for large values of Mi the two different calculations differ 



by many orders of magnitude. Using the simplest approximation, eqs |26| - |27| , 
it seems that larger values of Mi lead to larger values of Yb and that, in 
principle, one can get Yb as large as one desires. One can trace back this 
scaling to eq. ^ in section 3.2 and to the fact that eqs |2B| - P7| do not depend 
on the value of Mi. 



On the other hand, using eqs P^-P^l which take into account the sup 



presion of Yb for large values of Mi and rhi one gets a completely different 
picture. For large values of Mi, Yb is suppressed to negligible values and 
going to smaller values of Mi increases Yb- Note, however, that for the 
smallest values of Mi of the order of 0(10'^) [GeV] Yb stops growing and 
never reaches the experimentally preferred range of Yb ~ (0.5 — 1) x 10^^°. 
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Figure 8: Plots of Yg following eqs ^B|-^ (to the left) or according to eqs 
30| - |32| (to the right) for the different variants of model FC3. Line style as in 
fig-i 



All these variants of FC3 therefore fail the leptogenesis test. 

Since the simple approximation, eqs [26| - p7| pT[| , employed similarly by a 
number of authors |^, fails to take into account any Mi dependence of the 
dilution function its use would lead to the opposite conclusion. A careful 
treatment of d seems to be absolutely necessary for a reliable calculation of 
Yb and we stress again that also our treatment is still only approximate. 

In fig. H we plot Yb as defined in eqs ^ and ^ with ei estimated from eq. 
|5^ and with d calculated by a) the simple approximation, defined in eqs p6| - |27 
and b) our fit to the exact solution of the Boltzman equations [^, defined 
in eqs pil| - |52| . For both calculations we fixed rhi to rfii = 0.05 eV. For small 
values of Mi both approximations agree quite well, whereas for Mi larger 
than Ml ~ 10^ GeV the expectations from the different approximations 
differ by many orders of magnitude. 

One can understand the failure of the models FC3 to produce the correct 
value of Yb on the basis of the discussion presented in section 2.3. Q From 
the analytic estimates presented in section 2.3 one finds that in models where 
the dominant right-handed neutrino is the lightest, rhi depends on the same 
combination of Yukawas as the value of the heaviest neutrino mass, fixed by 
the atmospheric neutrino mass scale. Thus, for these models rhi ~ ""^s ~ 0.05 
eV. At such large values of rhi, however, the dilution function d, see Fig. |], 

®For the analytic estimations of ei and rhi we assumed that the right-handed neutrino 
mass matrix is diagonal, which is approximately true for the variants of FC3. 
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Figure 9: Plots of Yb following eqs p6| - p7| (broken line) or according to eqs 
30| - |32| (full line) assuming rhi = 0.05 eV. For this plot we have assumed that 
ei scales as given by eq. E^. 



is heavily suppressed for values of Mi larger than Mi ~ 10^ GeV. Thus, 
although larger values of Mi lead to larger values of ei, see eq. 0, the 
price one hast to pay for such large masses in the dilution function always 
overcompensates and Yb in these models can never be larger than Yb ~ 10~^^ 
as is demonstrated in fig. ^. 

Since all models where the dominant right-handed neutrino is the lightest 
share the feature rhi ~ ms, we conclude that upon use of eqs pPfP^ they all 
fail the leptogenesis test. In the next section we will therefore study models 
in which the dominant right-handed neutrino is the heaviest. 



7 Leptogenesis prefers models where the dom- 
inant right-handed neutrino is the heaviest 

In the previous section we have seen that although leptogenesis is decoupled 
from the low energy neutrino observables, nevertheless leptogenesis is capable 
of resolving the ambiguities between classes of models which would otherwise 
lead to the same experimental predictions. As an example of the power 
of leptogenesis to give information about the high energy theory, in this 
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section we show that leptogenesis prefers models where the dominant right- 
handed neutrino is the heaviest one and discuss the imphcations of this for 
unified models. According to our analytic estimates we expect this class of 
models to yield a lightest right-handed neutrino mass which is lighter than 
in the previous case, and hence more acceptable from the point of view of 
the gravitino constraint. In addition these models may be more consistent 
with GUTs. 

As a first example of a case (b) model we consider the charge vector 

(/i, k, h, ni,n2, rig, a) = (-3, 1, 1, 9, 1, -1, 2), (68) 

which defines a new model called FC9. The charges in Eq]6^ lead to 



y;^^ ~ A^o A2 1 (69) 
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Now if we take A = vA, this leads to the contributions from the heaviest 
(dominant), intermediate, and lightest right-handed neutrino, respectively, 
to the effective Majorana matrix of the order of, 

/ A2 A2 A2 \ / A2 A^ A^ \ 

A2 A2 A2 + O A^ A^o Aio . (71) 
V A2 A2 AV V A^ A^o A^o / 

By inspection we see that the model predicts 6^12 ~ 1, 9is ~ A^, and, from 
the order A^ accuracy of the SRHND condition, R ~ A^. It may therefore be 
suitable for one of the large mixing angle solar solutions, either LMA or LOW. 
Assuming A = a/A, the lightest right-handed neutrino mass is predicted to 
be X' ~ A^°F, or X' ~ S.IO^^F ~ 10^ GeV, which is rather small. In order 
to increase X' we need to increase A. 




As seen from Fig. |T0| one can adjust the hierarchy in the right-handed 



sector by a rather small change in A. Going from A = vA ~ 0.47 to A = 0.55 
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Figure 10: Plots of (from top left to bottom right: SAtm, s©, R, Sc, (S) and 
Ml [GeV] for the model defined in eq. ^ Dashed line: A = V^, full line: 
A = 0.55, dash-dotted: A = 0.60. 



31 




10"' 10"^ 10"= 10"^ 10-'^ 10" 




10-15 10-^3 10"" 10"^ 

Yb 



Figure 11: Plots of e (to the left) and Yb (to the right) for the model defined 



in eq. p8l Line styles are as in fig. 110. 



(0.60) changes Mi from Mi ~ (few) 10^ GeV to Mi ~ 10^° (10") GeV. This 
way it is possible to achieve larger values of e and a value of Yb marginally 



consistent with experimental data as shown in Fig. |Tl|. Note, however, that 
the peaks in Yb for this model are still too small and the models survive the 
leptogenesis test only in the tails of the distributions. Although in principle 
in models of this kind (case (b)) there is no association of the parameter m{ 
in Eq.^ with a physical neutrino mass, in this case when we calculate fhi 
we must first rotate to the basis in which the right-handed neutrino mass 
matrix is diagonal. This will lead to larger values of rfji ~ 10~^ — 10~^ eV 
than would naively be estimated from Eq.|5I| 

This change in A also influences (although only rather weakly) the pre- 
ferred values of R. As shown in Fig. |10] this model tends to prefer values 
of R consistent with the LOW solution of the solar neutrino problem. As 
mentioned previously, the LOW solution really only makes sense within the 
framework of SRHND because of the large hierarchies of neutrino masses 
which would otherwise appear rather unnatural. 

One of the advantages of having the dominant right-handed neutrino 
as the heaviest is that leptogenesis may be achieved consistent with e ~ 
1, which allows the third (dominant, and heaviest) right-handed neutrino 
to be associated with the third family in unified models. An example of 
such a model was recently presented in the framework of a string-inspired 
SUSY Pati-Salam (PS) model [^. The model in [^ will not be repeated 
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Figure 12: Plots of (from top left to bottom right: SAtm, s©, R, Sc, (S) and 



Ml [GeV] for the Pati-Salam model discussed in eqs [72| - |73| . Note that this 
model is consistent with the LA-MSW solution to the solar neutrino problem. 
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here, but we would emphasise that it was deduced from an analysis of the 
quark and lepton masses and mixing angles without any consideration of 
leptogenesis, and therefore we find it somewhat remarkable that it leads to 



a baryon asymmetry of the correct order of magnitude. The model in [25 



leads to the following structure for the Yukawa and right-handed neutrino 
mass matrix: 



(72) 





/ A^-s 


^3.5 


^1.5 


V 


^6.5 


^3.5 


1 




1^^6.5 


^3.5 


1 


( A' 


A^ 


A^ 


M^|~ A^ 


A^ 


^3.5 




A^ 


^3.5 


1 



(73) 



The effective light Majorana matrix then has contributions from the third, 
second and first right-handed neutrinos of: 



m 



PS 
LL 



A3 Ai-5 Ai-5 
A^-^ 1 1 
A^-^ 1 1 



O 



A2 A2 A2 
A2 A2 A2 
A^ A2 A2 



-O 



( A^ A^ A^ 

A^ A^ A^ 

V A^ A^ A^ 



• (74) 



From the analytic estimates in 3.1 we expect this model to be consistent 
with the LMA MSW solution. This is explicitly demonstrated in Fig. |12. 



From the analytic estimates in 3.2 we also expect this model to give successful 



leptogenesis, and this is demonstrated in Fig. |T3| From the matrices given 
above and from the analytical estimates of eq, 
3/(167r) X A^ ~ 7 X 10~®, which within a factor of 
numerical calculation of ei in fig 



53| one expects that ei ~ 
2 or so agrees with the 



13. Note, that the resulting values of Yb 



in fig. |T^ are also consistent with the arguments in 3.2. In particular in such 
case (b) models where the dominant right-handed neutrino is the heaviest it 



is easier to avoid the gravitino constraint |^, although the values of Mi are 
still a bit on the high side, as seen in fig. W\ As in the previous model, in the 
diagonal right-handed neutrino basis we obtain a larger value of m^ ~ 10~^ 
eV than would naively be estimated from Eq.|72|. 

Note that supersymmetric models of case (b) have the feature that there 
is an order unity Yukawa coupling in the 23 position of the Yukawa matrix 
which leads to a large off-diagonal entry in the slepton mass matrix. This 
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Figure 13: Plots of e (to the left) and Yb (to the right) for the Pati-Salam 
model of [25]. 



leads to the striking signature of the lepton flavour violating (LFV) process 



T —>■ ^•y close to the experimental upper limit, as first pointed out in [^ 



In general LFV constraints provide an additional window into the see-saw 



matrices in supersymmetric models |27, 28 . 



8 Summary and Conclusions 

This paper represents the first study of leptogenesis based on hierarchical 
models of neutrino masses in which SRHND is used to generate the neutrino 
mass hierarchy. Such models have been shown to accomodate the presently 
favoured large solar angle solutions such as LMA and LOW []TU|, and in the 
case of the LOW solution where the neutrino mass hierarchy is large it would 
seem that SRHND is almost inevitable. So we would argue that, far from this 
analysis being restricted to a particular small class of models, it is in fact quite 
generally applicable to large classes of models in which the neutrino mass 
hierarchy is generated in a natural way without any fine-tuning. Therefore 
the above results should be regarded as being quite generally applicable to 
see-saw models containing a neutrino mass hierarchy. 

In presenting our analytic and numerical results we make a clear distinc- 
tion between the theoretically clean asymmetry parameter ei and the baryon 
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asymmetry Yb, for which we present and use a fit based on a Boltzmann 
analysis. We have presented analytic expressions for both the MNS parame- 
ters, extending the previously presented analytic results [^ to the complex 
domain, and for leptogenesis asymmetry parameter ei in the cases where the 
dominant right-handed neutrino is either the heaviest or the lightest. We 
have compared the analytic estimates to full numerical results for models 
based on U{1) family symmetry, and have performed a numerical scan over 
the unknown coefficients, and have seen that the peaks of the distributions 
in ei are in good agreement with the analytic results. Using the analytic and 
numerical approaches we then discussed leptogenesis decoupling and lepto- 
genesis discrimination. 

We have shown that quite generally there is a decoupling between the low 
energy neutrino observables and the leptogenesis predictions for ei. Thus lep- 
togenesis has nothing to tell us about which solar solution we would expect, 
and for example the LMA and the LOW solutions are equally acceptable, as 
indeed would have been the SMA solution were it not disfavoured by SNO 
and Super-Kamiokande. Furthermore the leptogenesis phase is independent 
of the measurable MNS phase, although the analytic estimates make it clear 
that since the two phases originate from the same Yukawa matrix, and even 
in some cases involve the phases of the same Yukawa couplings, the general 
expectation is that, barring cancellations, both sorts of phases should be of 
roughly the same order of magnitude. 

In going from ei to Yb one needs to make some assumptions concerning 
the cosmological history of the universe. In this paper we have assumed a 
standard hot big bang universe, which is equivalent to assuming a very high 
reheat temperature after inflation which is larger than the right-handed neu- 
trino masses. Within this standard cosmology the right-handed neutrinos 
are produced via their couplings to the thermal bath, yet they are required 
to decay out-of-equilibrium, leading to a rather narrow range of couplings 
rhi of the lightest right-handed neutrino consistent with successful leptogen- 
esis. For the calculation of Yb a correct treatment of the Boltzman equations 
describing the number evolution of the heavy right-handed neutrinos in the 



early universe is essential [IS]. We therefore have devised an empirical flt 



formula and compared it to the exact results [|19[ as well as to the simpler 



approximation pO, |2TI|. Although for a small range in mi and small values 



of the lightest right-handed neutrino mass the simple approximation pO| , ^ 



agrees reasonably with the exact result [|T^, for most parts of the param- 
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eter space it fails badly. Only if one takes into account the suppresion of 
the dilution factor d for larger values of fhi and Mi, either by solving the 
Boltzman equations numerically or by the use of our approximate fit for- 
mula, does one find reliable results. Without taking this effect into account 
we would have wrongly concluded that Yb can get as large as Yb ~ 10~^ in 
some models, whereas with our more refined treatment we find that Yb is 
always Yb < 10~^^, if the dominant right-handed neutrino is the lightest. 

Based on the above analysis of Yb we have shown that leptogenesis ex- 
cludes a large class of models where the dominant right-handed neutrino is 
the lightest one. The power of leptogenesis to resolve ambiguities between 
models which would otherwise lead to the same neutrino observables provides 
a welcome constraint on high energy theories. We have shown that models 
where the dominant right-handed neutrino is the heaviest are marginally con- 
sistent with the gravitino constraint and have studied an explicit example of 
a unified model of this type. We find it encouraging that a model which 
was written down to describe the fermion mass spectrum [^, including the 
neutrino spectrum and the LMA MSW solution, should be precisely of this 
kind and gives successful leptogenesis, subject to the uncertainties of our 
estimates discussed in section 2.3. 

Finally we should emphasise that our conclusions are based on the as- 
sumed cosmological history being the standard hot big bang with a high 
reheat temperature. One plausible alternative is to suppose that the reheat 
temperature is below 10^ GeV, but that heavier right-handed (s)neutrinos 
can be produced in sufficient numbers by preheating at the end of inflation 



2^. The preheating must efficiently produce right-handed (s)neutrinos with- 
out over-producing gravitinos, and this will depend on the precise details of 
the inflation model. A model of leptogenesis with a low reheat temperature, 
based on preheating of heavy right-handed sneutrinos, which does not suffer 



from the gravitino problem has been recently studied in detail in |^6[. The 
same Pati-Salam model has also been studied in this context [^ and inter- 
estingly the results for Yg are also consistent, within the large uncertainty, 
with the estimates given here, based on an entirely different cosmological 
history of the universe. 
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